TITLES 

TYPE: 

CAPACITY? 

TEMPORARY STOEAGE: 

ACCURACY: 

NUMBER OF WORDS: 

DESCRIFTION: 



PROGRAM PARAMETER: 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER LABORATORY 

RY 






Symmetric Matrix Inversion 

Closed with 3 preset parameters 

2 < n < 136 where n ia the order of the matrix 

- 9 

A fimction of the order and conditioning of the matrix 

786 words on drum 

226 + 5n words in Williams Memory 

The routine inverts a syjtoetric matrix on the entry 



q 


TI nF 




50 q 


q+l 


26 (IN) 




00 X + pF 


q+2 


S F 




-, » 



There are several input options and output options. The 
input matrix nay be square or lower triangular and may be 
read in from the drum or from tape. The output may be punched 
or stored on the drum. If stored on the drim^ it will be 
in the form of a lower triangular matrix. If punched^ it 
may be in the form of a square^ lower triangular or a column 
of the diagonal entries. The input matrix can be scaled 
by 10*^ and the punched output can have the number of decimal 
places designated. 

1) T designates type, of matrix input 
S - triangular 

5 ^ square 

2) I designates input option 

= tape input (see note l) 

9 " drum input from drum location specified by S5 
5) n is the order of the matrix 
\') X is the output option 

t-i punch a lower triangular matrix 

1024 - punch a square matrix 

20^8 " pynch a column of the diagonal entries 



= 2 - 



EXAMPLES: 



THE INPUT MTEIX: 



5072 - store lower triangular matrix on drum beginning 
at location specified by S5. 

5) p is the number of decimal places in the punched output 

6) S is the scaling parameter for the input matrix 
y - positive scaler (see Note 2) 

LL 4096 - [y[ " negative scaler 
1) Input a square matrix of order 20 from the drum and punch 

the output as a triangular matrix to 6 places. Scale 

2 
the input by 10 o 

The entry is 



q. 


59 20F 




50 q 


q+l 


26 (in) 




00 6f 


q+2 


00 2F 







2) Input a triangular matrix of order kO from tape. Scale 

the input by 10 and punch the inverse as a square matrix 
to 5 places. 

The entry is 



q. 



q+1 



q+2 



SO 40F 
50 q 
2k (in) 
00 102 9F 



LL i|-095F 



The input matrix is in the form of signed fractions. It may 
be read in from the drum beginning at location S5 or it may 
be read in from tape. 

TAPE IHPUlP 

1) Triangular case. The matrix is punched by rows and termin= 
ated by a J, 

2) Square case. The matrix is punched by rows with an N 
terminating each row. Terminate the matrix with a Jo . 



- 



THE OUTPUT MA.TRIX; 



PRESET PARAMETERS: 



PROGRAMMED STOPS; 



DURATION; 



NOTE 1; 



The output is punched as signed fractions. 

TAPE OUTPUT 

1) Lower triangular caseo The matrix is punched hy rows. 
2 carriage returns and line feeds separate each row» 
The matrix is terminated by a Jo 

2) Square case» The matrix is punched by rows with an N 
terminating each rowo The matrix is terminated by a J. 

3) Diagonal caseo A single column of diagonal elements 
are punched o 

A p-place number follows the output in the above cases » 

This designates the location of the decimal point by locating, 

the point to the right of the column in which a 1 appears , 

v+1 
The output must be rescaled by 10'' 

The routine also can store the output as a lower triangular 

matrix on the drum beginning at location S5» 

Before this routine is read in^ the following parameters 

should be stored in location 3} ^3 and 5- They are of the 

form OOF OOvF 

V / S5 - location of this routine in Williams Memory 

Sk - " " 'a".:: " ." drum 

S5 - " " matrix in drum 
There are 5 programmed stops * Moving the white switch up and 
down at these points will allow the program to exit from 
this routine o 

FFOOO - lyfe trix is singular 
FFOOl - " " over-scaled 
FF002 - Diagonal element > 1 = 2°^^ 

The computation time is [ oOO^n-"^ + o08n + .In] seconds where 
n is the order of the matrix. 

The left hand order of q+1 can be changed to 2^ (IN) during 
tape input == The matrix may then be read in with a black 
switch start 6 



- k - 



NOTE 2 



The input matrix shoiild be scaled by 10^ s.Q that the largest 
elementLis < 1/2 in machine representation^ Y is any integer, 
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ORDER 
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00 226 








S3 02K 




226 


(Y) 


JO 18 


1 






S3 5OL 


y store part 1 on drum 


227 




26 (Yl) 








00 s^ 




228 




00 13^F 








26 999F 




229 




JO S3 
50 3L 


-V 


230 




26 (Yl) 








00 135S^ 


store part 2 on drum 


251 




00 151F 








26 999F 


y 


232 




JO S3 
50 6l 


^ 


233 




26 (Yl) 








.00 287S^ 


Store part 5 on drum 


23^ 




00 153F 








26 999F 


J 


235 




JO S3 
50 9L 




236 




26 (Yl) 








00 ii■i^-lS4 


Store part h on drum 


237 




00 88f 
26 999F 


. 


238 




JO 88S3 
50 12L 




239 




26 (Yl) 








00 530S4 


Store part 5 on drum 


2k9 




00 26F 
26 999F 


/ 


2i^-l 




JO 88S3 
50 15L 


■\ 


2^2 




26 (Yl) 


. 






00 557Si^ 


Store part 6 on drum 



LOCATION 




ORDER 
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2k3 




00 22F 




2kk 




26 999F 
JO 88S5 
50 18L 


J 


2^5 




26 (Yl) 


y 






00 58os4 


store part 7 on drum 


2k6 




00 lOF 








26 999F 


J 


2k^ 




JO S3 
50 21L 


\ 


248 




26 (Yl) 








00 591S4 


*■ Store part 8 on drum 


2^9 




00 56F 








26 999F / I 


250 




J-O S3 
50 2ij-L 


-\ 


251 




26 (Yl) 








00 6i^8si^- 


Store part 9 on drum 


252 




00 137F 








26 999F >^ 






00 153S3 








02 K 




153 


(N) 


00 F 








00 F 


Constants 


l^k 


(P) 


00 F 
00 F 




155 


(X) 


00 F 
00 F 




156 


(s) 


OOjF 
00 F 




157 


(I) 


00 F 
00 F 




158 


(MS) 


00 F 
00 F 




159 


(z) 


00 F 








00 F 





LOCATIOl^ 


ORDER 
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160 


(VI) 


7J 226S3 
i^o 226S3 




161 


(V2) 


00 F 
00 F 




162 


(V3) 


Ih F 
ho F 




165 


(1) 


00 F 
00 IF 




iGk 


(1-1) 


00 IF 
00 IF 




165 


(+) 


00 F 
00 lOF 




166 


(") 


00 F 

00 1000 0000 0000 J 




167 


(+1) 


00 F 

00 9999 9999 9999 J 




168 


(3^) 


00 F 
00 3^F 




169 


(CI) 


00 F 
00 F 




170 


(Dl) 


26 139S3 
00 F 

00 K 




171 


(mtk) 


50 18S3 
50 L 




172 




26 (Yl) 








00 S^ 


' Playback part 1 


173 




00 13-^F 

26 18S3 




17^ 




50 S3 
50 3L 




175 




26 (Yl) 


- 






00 135S4 


1 Playback part 2 


176 




00 151F 


) 






26 S3 
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177 




50 S3 


"1 






50 6l 




178 




26 (Yl) 








00 287S4 


f Playback part 3 


179 




00 153F 
26 S3 


J 


180 




50 S3 
50 9L 


\ 


181 




26 (Yl) 


V. 






^^ 1.1. to'' 


/ Playback part 4 


182 




00 88F 
26 S3 


J 


183 




50 S3 
50 12L 


] 


184 




26 (Yl) 








00 648S4 


Playback part 9 


185 




00 137F 




186 


(n) 


22 F 
00 K 

00 18S3 
02 K 


/ 


18 




L5 (V2) 








42 43 (lA) 


Initialize addresses for (lA) 


19 




42 26 (lA) 
42 4l(lA) 




20 




00 20F 
46 6(1A) 




21 




F5 (V2) 
42 23 (lA) 




22 




42 24 (lA) 
42 12 (lA) 




23 




L5 (V3) 
42 6(1A) 




24 




F5 (V3) 
42 10 (lA) 




25 




42 78 (lA) 






- 


42 79 (lA) 
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00 K 




26 


(lA) 


L5 (1) 








40 2F 




27 




L4 (1) 

40 3F 




28 




LI (N) 
L4 (1) 




29 




4o (ci) 

40 7F 




30 




L5 (1) 

4o 6f 




31 




85 11^ 








00 S5 




32 




40 F 
L5 F 




35 




4o 4f 
22 8L 




34 




50 226S3 


"?V 






50 8L 


Playback one row of A 


35 




26 (Yl) 

00 2S5 




36 




00 2F 
LI F 




37 




l4 4f 
4o IF 




38 




L5 226S3 
4o F 




39 




L3 IF 
32 25L 




40 




L7 IF 
00 2py 




41 




46 20L 
46 24L 




42 




LI IF 
32 23L 




43 




l4 4f 








40 4f 
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kk 




L5 (Vg) 
i^-2 19L 




h^ 




k2 20L 
L5 F 




k6 




10 F 
40 F 




hi 




F5 20L 








LO 24L 




kQ 




36 25L 
F5 19L 




h9 




22 18L 








L5 F 




50 




10 F 
iK) F 




51 




i^l IF 
L5 226S3 




52 




i|0 F 








L7 F 




53 




L2 F 
32 30L 




5^ 




F5 IF 
i^O IF 




55 




L5 F 
10 IF 




56 




26 26l 
L3 IF 




57 




32 39L 
L7 IF 




58 




00 20F 
k6 36L 




59 




Fl 2F 
kO 5F 




6o 




L5 (VI) 
1+2 35L 




61 




4-2 36L 
L5 F 





.OCATIO Ill 
62 

6k 

65 
66 
61 
68 

69 
70 

71 
72 

73 
7^ 
75 
76 
77 
78 
79 



ORDER 
10 F 
kO F 
E5 5F 
i|-0 5F 
52 59L 
F5 55L 
22 3^L 
L7 ^F 
lA (1) 
00 20F 
k6 i^-5L 
L3 F 
36 124l 
50 226S3 

75 (+1) 

66 F 

S5 F 

i^-0 8f 

10 IF 

00 IF 

32 14-6L 

ko 226S3 

Lli 6f 

^0 5F 

F5 ^3L 

k2 5OL 

F5 ^2L 

1^2 5IL 

42 5ifL 

L5 F 

kO F 

L5 F 

40 9F 

36 53L 

50 8f 

79 F 



NOTES 
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KSL 5-14 



Test for singularity; Jump if yes 



J 



Form 



^ik 



11 



Form X • a 
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80 




lJ+ 9F 
ho F 


J 


81 




F5 5F 








i^O 5F 




82 




52 58L 








F5 50L 




83 




k2 50L 
F5 51L 




Sk 




22 49L 

Fl 2F 




85 




40 F 
L5 (VI) 




86 




42 6lL 
i|-2 62L 




87 




L7 5F 
L2 F 




88 




52 65L 
L7 F 




89 




i+0 5F 
.F5 F 




90 




40 F 
56 66l 




91 




F5 6IL 
26 6OL 




92 




Fl 2F 
ko 9F 




93 




LL 5F 








52 lllL 


Test if ^. . > 1/2; jump if no 


9^ 




F5 F 
k2 74l 




95 




FO (1) 

LO IF 




96 




40 F 
L5 F 




97, ' 




52 78L 








L5 (VI) 
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98 


42 75L I 






42 75L 




99 




50 (1) 

L5 F 




100 




10 2F 
00 IF 




101 




52 75L 
40 F 




102 




F5 9F 








52 78L 




105 




40 9F 








F5 75L 




104 




26 72L 
L5 F 




105 




L4 F 
40 F 




106 




JO 226S5 

50 8OL 


■s 


107 




26 (Yl) 
00 2S5 


store row of L 


108 




00 2F 




109 




F5 7F 
52 95L 
40 7F 


^ 


110 




L5 lOL 
L4 (1=1) 




111 




40 lOL 
46 82L 




112 




42 78L 
42 79L 




115 




F5 2F 
40 2F 




114 




L4 (1) 

l4 9L 




115 




4o 9L 
40 81L 
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116 




F5 12L 
k2 12L 




117 




i|2 23L 
k2 2iH 




118 




F5 6f 
ko 6f 




119 




22 8L 
F5 (CI) 




120 




36 115L 
F5 6l 




121 




k2 6l 

L^ (1) 




122 




k2 lOL 
k2 78L 




123 




k2 79L 
F5 (V2) 




12^ 




^2 12L 
k2 2^L 




125 




i|-2 2^L 

L5 3F 




126 




i|0 2F 
L^ (1) 




127 




ko 3F 
L^ (C2) 




128 




^0 (C2) 
^0 9L 




129 




i^-0 8IL 
L5 3F 




130 




00 20F 
k6 lOL 




131 




^6 82L 
F5 25L 




132 




^2 25L 
il-2 i^2L 




133 




k2 k6L 






, 


L5 3F 





LOCATION 



13^ 
135 
136 
137 



1 qQ 



139 

1^0 

1^3 
150 



ORDER 




F5 F 

M-U r 

L7 5F 
50 F 
00 IF 
^0 5F 
26 67L 

L5 (1) 
4o 3F 
LI (n) 
^0 F 
40 5F 
L5 (V3) 
^2 119L 
k2 120L 
L5 (MS) 
L^ F 
32 121L 
L7 F 
^0 (MS) 
F5 F 

32 3(mtr) 

^0 F 

F5 119L 

26 II8L 

FF F 

22 12(MTR) 
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FFOOO: matrix singular, -wh. sw. to 
continue. 
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151 


'(C2) 


26 (Yl) 
00 2S5 




152 




22 (y) 
26 IN 

00 K 

00 S302K 









L5 (V2) 








il-2 18(2A) 


Initialize addresses for (2A) 


1 




k2 lOil-(lT) 
k2 105 (IT) 




2 




L5 (V3) 
h2 6(2A) 




3 




i^-2 33(2A) 
^0 53(IT) 




ii- 




h2 9^(IT) 
i^2 95(IT) 




5 




L5 (N) 

^2 art IT) 




6 




00 20F 
k6 87(iT) 

00 K 




7 


(2A) 


F5 (MB) 
00 20F 




8 




k6 2L 
L5 (+1) 




9 




10 F 
00 IF 




10 




ho 6f 

L5 (VI) 




11 




^2 13L 
^2 1^4-1 




12 




L5 3F 
Li^ 13L 
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13 




kO k¥ 

L5 F 






Ik 




L^ (m) 

kO 2F 






15 




50 226s 3 


Playback ro-w of L 








50 8L 




16 




26 (Yl) 
00 S5 


^ 




17 




00 IF 
L3 2F 






18 




36 18L 
L7 2F 






19 




00 20F 

kS ikL 






20 




01 IF 
L5 F 






21 




10 F 
kO F 






22 




F5 13L 
k2 13L 






23 




k2 l^L 
LO kF 






2k 




36 18L 
22 13L 






25 




L5 226s 3 
^0 F 






26 




L5 6f 
kO 226s 3 






27 




JO 226s 3 

50 20L 


•• 




28 




26 (Yl) 

00 S5 


Store row of L (rescaled) 




29 




00 IF 

F5 5F 







LOi^AXXQI 
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30 




36 33L 

ko 5F 




31 




F5 6l 
k2 6l 




32 




F5 3F 
ko 3F 




33 




Li^- 9L 
^0 9L 




34 




40 21L 
L5 lOL 




35 




Lk (1-1) 
.i^6 lOL 




36 




A6 PPT, 

L5 i8l 




37 




Lk (1-1) 

ko i8l 




38 




F5 19L 
i^2 19L 




39 




22 3L 
50 F 




i^-0 




LI (MS) 
ko F 

00 K 




ki 


(IT) 


L5 (C2) 








4o i8l 


Form L routine 


k2 




ko 77L 

LI (N) 




h3 




^0 2F 
Lk (1) 




hk 




il-0 IF 
L5 87L 




^5 




LO (1-1) 








i^O 3F 


,^ 
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ke 




LO (1-1) 
Lk 22L 




hi 




40 PPT, 
LI IF 




k^ 




Lk 88l 
40 88l 




k9 




L5 70L 
LO 2F 




50 




4o 87L 
L5 3F 




51 




l4 53L 
40 53L 




52 




42 i6l 

Ll 2F 




53 




4o 3F 
4o 4f 




54 




00 20F 
46 19L 




55 




46 78L 

F5 (MS) 




56 




00 20F 

46 30L 




57 




Ll (MS) 
40 F 




58 




50 226S3 








50 17L 


1 PlaylDack row of L 


59 




26 (yi) 
00 F 




60 




00 F 

L5 i8l 




61 




LO 91L 

l4 90L 




62 




4o 5F 








4o 6f 
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63 




LI 226S3 

ko 226S5 


Form ^^_^^r^ 




Gk 




F5 IF 
52 76L. 






65 




L5 22L 
LO (1) 






G6 




k2 29L 
k2 3^L 






67 




k2 57L 
L5 86l 






68 




ko 57L 
ij-1 F 






69 




^1 7F 
in 8f 






70 




50 7F 
L5 F 






71 




10 F 
00 IF 






72 




26 55L 
ifO 5F 






75 




85 IIF 
00 F 






7^ 




ko 9F 
50 9F 






75 




L5 7F 

Ik F 






76 




^0 9F 
S5 F- 


Form £^j" 




77 




ko 7F 
L5 9F 






78 




10 IF 






79 




00 IF 
L^ -F 
^0 F 


ii 




80 




LL F 










36 Wi. 







LOCATION ' 
81 

82 

83 
Qk 

85 
86 

87 
88 

89 
90 
91 
92 
93 
9^ 
95 
96 
97 



ORDER 
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L5 F 
10 IF 
kO F 
F5 8F 

ko 8f 
L5 37L 
I^ (1-1) 
i|-6 37L 
F5 54l 
42 34L 
LO 88l 
32 6hL 
L5 5F 
ko 32L 
lA 3F 
ko 5F 
F5 3F 
22 31L 
L5 F 
00 IF 
ko F 
F5 8F 
36 57L 
26 66l 
L3 8f 
36 57L 
lA F 
40 F 
LI 8F 
50 (1) 
00 20F 
Lk 3OL 
46 3OL 
50 F 
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98 




LI F 
^0 F 




99 




L5 89L 
LO 57L 




100 




32 74l 
L5 6f 




101 




LO kF 
40 5F 




102 




ko 6f 

L5 ^F 




103 




LO (1) 
^0 3F 




lOl^ 




L5 57L 




105 




22 24L 
50 7F 




lo6 




LI 8F 
LO (MS) 




107 




ko 8f 

LL F 




108 




36 49L 
L7 8f 




109 




00 20F 

k6 71L 




110 




F5 29L 
k2 7OL 




111 




^2 7IL 
L5 226S3 




112 




10 F 
^0 F 




113 




F5 70L 

LO 87L 




lli^ 




36 52L 






< 


F5 70L 
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ORDER NOTES PAQS 19 KSL 5»1^ i 


115 




22 69L 
L5 22L 




116 




LO (1-1) 
kO 22L 




117 




JO 226S3 

50 76L 


1 


118 




26 (Yl) 

00 F 


Store row of L 


119 




00 F 








F5 IF J 


120 




36 92L 
kO IF 




121 




L5 18L 
lh 2F 




122 




ho 18L 
kO 77L 




123 




F5 2F 
kO 2F 




idb' 




L5 88l 
LO (1) 




125 




k2 88l 
LI (1-=1) 




126 




26 lOL 
50 F 




127 




10 IF 
00 IF 




128 




00 F 
00 F 




129 




L5 7F 
7^ 2i6S3 




130 




LI F 








ko 226S3 




131 




85 IIF 




- 




00 F 
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li^9 


k2 lOi^L 










k2 105L 






150 




22 10i|-L 
50 F 






151 




22 3(Y) 

26 IN 

00 K 

00 S502K 




\ 







L5 (V2) 










^2 4(d) 


Initialize addresses 


for (3A) 


1 




42 6(D) 
k2 22(d) 






2 




42 23(D) 
L5 (V3) 






3 




42 (5A) 
92 l47F 

00 K 






1^- 


(3A) 


L5 (MS) 
l4 F 






5 




L4 (1) 

42 7L 






6 




L5 (VI) 

42 6l 






7 




42 8L 
22 4l 






8 




50 226S3 
50 4l 


1 

/- Playback row of L 




9 




26 (Yl) 
00 S5 






10 




00 IF 
50 F 






11 




7J (^) 










10 F 
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12 




00 IF 
kO F 




13 




F5 8l 








k2 6l 




Ik 




42 8L 






LO 2i^L 




15 




32 12L 
22 6L 




16 




JO 226S3 
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